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Abstract 

An initial- value problem (IVP) for arbitrary small three-dimensional vorticity 
perturbations imposed on a free shear flow is considered. The viscous perturbation 
equations are then combined in terms of the vorticity and velocity, and are solved 
by means of a combined Laplace-Fourier transform in the plane normal to the 
basic flow. The perturbations can be uniform or damped along the mean flow 
direction. This treatment allows for a simplification of the governing equations such 
that it is possible to observe long transients, that can last hundreds time scales. 
This result would not be possible over an acceptable lapse of time by carrying 
out a direct numerical integration of the linearized Navier-Stokes equations. The 
exploration is done with respect to physical inputs as the angle of obliquity, the 
symmetry of the perturbation and the streamwise damping rate. The base flow 



is an intermediate section of the growing two-dimensional circular cylinder wake 
where the entrainment process is still active. Two Reynolds numbers of the order 
of the critical value for the onset of the first instability are considered. The early 
transient evolution offers very different scenarios for which we present a summary 
for particular cases. For example, for amplified perturbations, we have observed 
two kinds of transients, namely (1) a monotone amplification and (2) a sequence 
of growth - decrease - final growth. In the latter case, if the initial condition is 
an asymmetric oblique or longitudinal perturbation, the transient clearly shows an 
initial oscillatory time scale. That increases moving downstream, and is different 
from the asymptotic value. Two periodic temporal patterns are thus present in the 
system. Furthermore, the more a perturbation is longitudinally confined the more 
it is amplified in time. The long-term behavior of two-dimensional disturbances 
shows excellent agreement with a recent two-dimensional spatio-temporal multiscale 
modal analysis and with laboratory data concerning the frequency and wave length 
of the parallel vortex shedding in the cylinder wake. 
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1 Introduction 

Recent shear flows studies ([I], [2], [3]) have shown the importance of the early time dy- 
namics, that in principle can lead to non-linear growth long before an exponential mode 
is dominant. The recognition of the existence of an algebraic growth, due - among other 
reasons - to the non-orthogonality of the eigenfunctions (|3]) and a possible resonance be- 
tween Orr-Sommerfeld and Squire solutions ([5]), recently promoted many contributions 
directed to study the early-period dynamics. For fully bounded flows works by [3J, [UJ, 
[Z], [E], [9], [ID], and for partially bounded flows works by [TT] . [T2] . [13] . [H] can be cited. 
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As for free shear flows, the attention was first aimed in order to obtain closed-form so- 
lutions to the initial- value inviscid problem ([13], [IS]) and was successful by considering 
piecewise linear parallel basic flow profiles. 

An interesting aspect observed in the intermediate periods is that the maximal ampli- 
fication is generally associated to oblique disturbances, that, as a consequence, potentially 
can promote early transition, see e.g. [13] ■ In fact, these perturbations, which are asymp- 
totically stable at all Reynolds numbers, are the perturbations best exploiting the energy 
transient amplification. 

In this work we consider as a prototype for free shear flow the two-dimensional wake 
past a bluff body. The wake stability has been widely studied by means of modal anal- 
yses (e.g. [17], [H], [19], [20] ). However, in this way only the asymptotic fate can be 
determined, regardless of the transient behavior and the underlying physical cause of any 
instability. 

In this work we adopt the velocity-vorticity formulation to evaluate the general initial- 
value perturbative problem. This method was proposed by Criminale, Drazin and co- 
authors in the years 1990-2000 ([2], [H], [IB] , [BJ, [H]). In synthesis, the variables are 
Laplace-Fourier transformed in the plane normal to the basic flow. Afterward, the result- 
ing partial differential equations in time are integrated numerically. This procedure allows 
for completely arbitrary initial expansions by using a known set of functions (Schauder 
basis in the L 2 space) and yields the complete dynamics — the early time transients and 
the asymptotic behavior (up to many hundred time scales) — for any disturbance. The 
long term dynamics would not, in fact, have been easily recovered by using the alternative 
method of the direct numerical integration of the linear equations since the integration 
over a range of time larger than a few dozen basic time scales is not feasible. 

The base flow model that we employ includes the wake transversal velocity and thus 
the nonlinear and diffusive dynamics that are responsible for the growth of the flow and 
the associated mass entrainment. We consider the first two order terms of the analytical 



Navier-Stokes expansion obtained by Tordella and Belan (2002, 2003) ([21], [22]), see {2.1 
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In particular, we consider the longitudinal component of such an expansion solution, the 
problem is parameterized by Xq, the longitudinal coordinate, and the Reynolds number 
Re. 

We use a complex wavenumber for the disturbance component aligned with the flow 
so that longitudinal spatially damped waves are represented. It should be observed that 
a longitudinal spatial growth could not be considered physically admissible as an initial 
condition since the energy density of the initial perturbation would be infinite. In the 
context of the initial-value problem, this is an innovative feature adopted to introduce a 
possible spatial evolution (damping) of the perturbative wave in the longitudinal direc- 
tion. The perturbative equations are numerically integrated by the method of lines. The 



equations formulation and initial and boundary conditions are presented in £2.2 



The perturbation evolution is examined for the base flow configurations corresponding 
to the Reynolds numbers of 50, 100, and for a typical section, xq = 10, of the intermediate 
region of the flow where the entrainment process is active. A comparison with a base flow 
far field configuration, xq = 50, is also proposed. The normalized perturbation kinetic 
energy density is the physical quantity on which the transient growth is observed (see 



£3.1). To determine the temporal asymptotics of the disturbance an equivalent of the 



modal temporal growth rate is introduced (see £3.3). 

In the case of longitudinal disturbances, comparison with recent spatio-temporal multi- 
scale Orr-Sommerfeld analysis ([23], [21]) and with laboratory experimental data Q27J) 
is carried out, see §3.3[ As noted, the initial conditions posed are arbitrary. As far as 
the modal theory is concerned, the agreement is excellent both for the frequency, defined 
as the temporal derivative of the perturbation phase, and the temporal growth rate. It 
is also in quantitative agreement with respect to the laboratory data. The experiment 
shows couples of pulsation and wavelength of the cylinder vortex shedding that are close 
to those yielded by the IVP analysis when the wavenumber is that where the growth rate 
is maximum. Conclusions are given in §|4| 
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2 Initial-value problem 



2.1 Base flow 

The base flow is considered viscous and incompressible. To describe the two-dimensional 
growing wake flow, an expansion solution for the Navier-Stokes two-dimensional steady 
bluff body wake ([22], [21]) has been used. The x coordinate is parallel to the free stream 
velocity, the y coordinate is normal. This approximated analytical Navier-Stokes solution 
incorporates the effects due to the full non linear convection as well as the streamwise 
and transverse diffusion. The solution was obtained by matching an inner Navier-Stokes 
expansion in terms of the inverse of the longitudinal coordinate x (x~ n ^ 2 , n = 0, 1, 2, . . .) 
with an outer Navier-Stokes expansion in terms of the inverse of the distance from the 
body. 

Here we take the first two orders (n = 0, 1) of the inner longitudinal component of the 
velocity field as a first approximation of the primary flow. In the present formulation the 
near-parallel hypothesis for the base flow, at a longitudinal position x = xq, is made. The 
coordinate Xq plays the role of parameter of the steady system together with the Reynolds 
number. The analytical expression for the profile of the longitudinal component is 

_Rey2 

U(y; x , Re) = 1 - adx 1/2 e 4 x o (1) 

where a is related to the drag coefficient (a = - (Re / 'ix Y ' ' 2 cd (Re)) and C\ is an integration 
constant depending on the Reynolds number. As said in the introduction, this two terms 
representation is extracted from an analytical asymptotic expansion where the velocity 
vector and the pressure are determined to the forth order. It should be observed that 
the transversal velocity component V first appears at the third order (n = 2), while the 
pressure only at the forth order (n = 3). Up to the second order, the field is thus parallel. 
Beyond the second order the analytical expression becomes much more complex, special 
functions as the confluent hypergeometric functions plays a role ([21]) associated to the 
deviation from parallelism. By changing the Xq values, the base flow profile ([!]) will locally 
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Figure 1: Wake schematic. Profile Uf(y; xq, Re) in the intermediate (xq = 10) and far 
(xq = 50) wake for different Reynolds numbers, Uf is the free stream velocity. The 
diameter of the cylinder is out of scale (three times) with respect to the wake profiles. 

approximate the behavior of the actual wake generated by the body. Here, the region 
considered, if not otherwise specified, is fixed to a typical section, xq = 10D (where D is 
the spatial scale of the wake) of the intermediate wake. The term intermediate is used 
in the general sense as used by [2B]: 'intermediate asymptotics are self-similar or near- 
similar solutions of general problems, valid for times, and distances from boundaries, 
large enough for the influence of the fine details of the initial/or boundary conditions 
to disappear, but small enough that the system is far from the ultimate equilibrium 
state...'. The distance beyond which the intermediate region is assumed to begin varies 
from eight to four diameters D for Re e [20,40] (see [21], [22]). Base flow configurations 
corresponding to a Re of 50, 100 are considered. In Figure [T] a representation of the wake 
profile at differing longitudinal stations is shown. 
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2.2 Formulation 

By exciting the base flow with small arbitrary three-dimensional perturbations, the con- 
tinuity and Navier-Stokes equations that describe the perturbed system are 

du dv dw 
dx dy dz 

du ^du „dU dp 1 
dt dx dy dx Re 

®v y.dv ^dp _ 1 

dt dx dy Re 

dw ^ ^dw ^ dp _ 1 y2 ~ /g\ 
dt dx dz Re 

where (u(x,y, z,t), v(x,y, z,t), w(x,y, z,t)) and p(x,y,z,t) are the perturbation velocity 
components and pressure respectively. The independent spatial variables z and y are 
defined from — oo to +oo, while x is defined in the semispace occupied by the wake, from 
to +oo. All physical quantities are normalized with respect to the free stream velocity 
Uf, the body scale D and the density. By combining equations ^ to ^ to eliminate 
the pressure, the linearized equations describing the perturbation dynamics become 

,d TT d ,_ 9 _ dv d 2 U 1 _ 4 „ . . 

,d TT d dv dU 1 _ 9 „ ._. 

where u y is the transversal component of the perturbation vorticity field. By introducing 
the quantity f , that is defined by 

V 2 v = f (8) 

we obtain three coupled equations ^ and Equations ^ and ^ are the Orr- 
Sommerfeld and Squire equations respectively, from the classical linear stability analysis 
for three-dimensional disturbances. From kinematics, the relation 

~ du z doj x 
dx dz 
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physically links together the perturbation vorticity components in the x and z directions 
(uj x and ui z respectively) and the perturbed velocity field. By combining equations (|6| 
and g then 

dT TT dT dvd 2 U 1 2 ~ 

— \-u = — v 2 r (10) 

dt dx dx dy 2 Re 
which, together with (J7|) and (|8p, fully describes the perturbed system in terms of vor- 
ticity. This formulation is a classical one. Alternative classical formulations, as the 
velocity-pressure one, are in common use. We chose this formulation because the vortic- 
ity transport and diffusion is the principal phenomenology for the dynamics of a wake 
system. For piecewise linear profiles for U analytical solutions can be found. For con- 
tinuous profiles, the governing perturbative equations cannot be analytically solved in 
general, but may assume a reduced form in the free shear case Q26J). 
Moreover, from the equations (J7]), ^ and (10), it is clear that the interaction of the 



mean vorticity in ^-direction (Q z = —dll/dy) with the perturbation strain rates in x and 
z directions (|| and |^ respectively) proves to be a major source of any perturbation 
vorticity production. 

The perturbation quantities are Laplace and Fourier decomposed in the x and z 
directions, respectively. A complex wavenumber a = a r + ioti along the x coordinate 
as well as a real wavenumber 7 along the z coordinate are introduced. In order to have 
a finite perturbation kinetic energy, the imaginary part CKj of the complex longitudinal 
wavenumber can only assume non-negative values. In so doing, we allow for perturbative 
waves that can spatially decay (a, > 0) or remain constant in amplitude (oj = 0). The 
perturbation quantities (v, T, u y ) involved in the system dynamics are now indicated as 
(v, f , u y ), where 

/+00 r+00 
/ g(x } y,z } t)e- iax - ilz dxdz (11) 
-00 Jo 

indicates the Laplace-Fourier transform of a general dependent variable in the a — 7 
phase space and in the remaining independent variables y and t. The governing partial 
differential equations are 
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d 2 v 

-4 - (k 2 - ol + 2ikcos((p)a i )v = f (12) 
ay 2 



dT - d 2 U 

— = — (ikcos((j)) — a^UT + (ikcos(4>) — ai)——v 
at ay 2 

1 d 2 f 

+ -^[^-{k 2 -a 2 + 2tkcos{<P)a i )r] (13) 
He oy z 

= — (ikcosicj)) — a,i)UCj v — iksin(<f)) < ^-v 
at ay 

1 d 2 Cj. 

R~e [ ~dy 2 



+ —\^-f~{k 2 -a 2 + 2ikcos{(t))a l )Cj y \ (14) 



where = tan 1 ( r y /a r ) is the angle of obliquity with respect to the x-y physical plane, k 



a 2 + 7 2 is the polar wavenumber and a r = kcos((f)), 7 = ksin((J)) are the wavenumbers 
in x and z directions respectively. The imaginary part «j of the complex longitudinal 
wavenumber represents the spatial damping rate in the streamwise direction. In figure [2] 
the three-dimensional perturbative geometry scheme is depicted. 



From equations ( 12 )-( 14 ), it can be noted that there can neither be advection nor pro- 
duction of vorticity in the lateral free stream. The vorticity can only be diffused since 
only the diffusive terms remains in the limit when y — > 00. Perturbation vorticity van- 
ishes in the free stream, regardless if it is initially inserted there (if inserted, vorticity is 
finally dissipated in time when y — > 00). This means that the velocity field is harmonic 
as y — > 00. 



Governing equations (12), (13) and (14) need proper initial and boundary conditions 
to be solved. Among all solutions, those whose perturbation velocity field is bounded in 
the free stream are sought. Periodic initial conditions for 

r = t!-!t - (k 2 - a 2 + 2ikcos((j))a i )v (15) 
oy 2 

can be cast in terms of a set of functions in the L 2 Hilbert space, as 

v(0, y) = e~ iy ~ yo)2 cos(n (y - y )) or v(0, y) = e~ {y ' yo)2 sm(n (y - y )), 
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Figure 2: Perturbation geometry scheme. 

for the symmetric and the asymmetric perturbations, respectively. Parameter no is an 
oscillatory parameter for the shape function, while y$ is a parameter which controls 
the distribution of the perturbation along y (by moving away, or bringing nearer, the 
perturbation maxima from the axis of the wake). The trigonometrical system is a 
Schauder basis in each space L p [0, 1], for 1 < p < oo. More specifically, the system 
(1, sin(n y), cos(n y), . . .), where n = 1,2,..., is the Schauder basis for the space of 
square-integrable periodic functions with period 2ir. This means that any element of the 
space L 2 , where the dependent variables are defined, can be written as an infinite linear 
combination of the elements of the basis. 

The transversal vorticity u y is chosen initially equal to zero throughout the y domain in 
order to ascertain which is the net contribution of three-dimensionality on the transver- 
sal vorticity generation and temporal evolution. However, it can be demonstrated that 
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the eventual introduction of an initial transversal vorticity does not actually affect the 
perturbation temporal evolution. 

Once initial and boundary conditions are properly set, the partial differential equa- 



tions (12)-(14) are numerically solved by the method of lines. The spatial derivatives 
are centre differenced and the resulting system is then integrated in time by an adap- 
tative multi-step method (variable order Adams-Bashforth-Moulton PECE solver). The 
transversal computational domain is large thirty body scales. By enlarging the com- 
putational domain to 50 and 100 body scales the results vary on the third and fourth 
significant digit, respectively. 

2.3 Measure of the growth 

One of the salient aspects of the initial-value problem is to observe the early transient 
evolution of various initial conditions. To this end, a measure of the perturbation growth 
can be defined through the disturbance kinetic energy density in the plane (a, 7) (see e.g 

!) 



1 1 f+V* 2 2 2 

e(t;a,j,Re) = -— / (|u| + \v\ +\w\)dy 
* ^Vd J -yd 

= l^r^^r\\^\ 2 + \^ + l 2 \\v\ 2 + H\ 2 )dy, (16) 

2 2y d \a 2 + 7 2 \J-y d dy 

where 2yd is the extension of the spatial numerical domain. The value ya is defined so that 
the numerical solutions are insensitive to further extensions of the computational domain 
size. Here, we take = 15. The total kinetic energy can be obtained by integrating 
the energy density over all a r and 7. The amplification factor G(t) can be introduced in 
terms of the normalized energy density 

G{t; a, 7 = — — r. 17 

e[t = 0; a, 7) 

This quantity can effectively measure the growth of a disturbance of wavenumbers (a, 7) 
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at the time t, for a given initial condition at t = (Criminale et al. 1997; Lasseigne e£ 
a/. 1999). 

The temporal growth rate on the kinetic energy r 

,1';^) = ™, *>o (is) 

is introduced in order to evaluate both the early transient as well as the asymptotic 
behavior of the perturbations (here it is the first moment of the perturbation which is 
assumed to asymptotically approach an exponential growth). Computations to evaluate 
the long time asymptotics are made by integrating the equations forward in time beyond 
the transient ([5J, [TTJ) until the temporal growth rate r asymptotes to a constant value 
(dr/dt < e, where e is of the order 10~ 4 ). The angular frequency (pulsation) u of the 
perturbation can be introduced by defining a local, in space and time, time phase tp of 
the complex wave at a fixed transversal station (for example y = 1) as 

v(y, t; a, 7 , Re) = A t (y; a, 7 , Re)e^ (19) 

and then computing the time derivative of the phase perturbation (p (|29j) 



Since (p is defined as the phase variation in time of the perturbative wave, it is reasonable 
to expect constant values of frequency, once the asymptotic state is reached. 



3 Results 

We present a summary of the most significant transient behavior and asymptotic fate of 
the three-dimensional perturbations. The temporal evolution is observed in the interme- 
diate asymptotic region of the wake, which is the region where the spatial evolution is 
predominant. It can be demonstrated that changing the number of oscillations no and 
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the parameter yo that controls the perturbation distribution along the y direction can 
only extend or shorten the duration of the transient, while the ultimate state is not al- 
tered. More specifically, if the perturbation oscillates rapidly or is concentrated mainly 
outside the shear region of the basic flow, for a stable configuration, the final damping is 
accelerated while, for an unstable configuration, the asymptotic growth is delayed. Thus, 
these two parameters are not crucial, because their influence can be recognized a priori. 
Therefore, in the following we use the two reference values, no = 1 and yo = 0, and focus 
the attention mainly on parameters such as the obliquity, the symmetry, the value of the 
polar wavenumber and the spatial damping rate of the disturbance. In particular, the 
polar wavenumber k changes in a range of values reaching at maximum the order of mag- 
nitude 0(1), according to what is suggested by recent modal analyses ([24], [23]). The 
order of magnitude of the spatial damping rate a, t varies around the polar wavenumber 
value. 

3.1 Exploratory analysis of the transient dynamics 

Figure [3] takes into account the influence, on the early time behavior, of the perturbation 
symmetry and of the wake region considered in the analysis, which is represented by 
the parameter xq. All the configurations considered are asymptotically amplified, but 
the transients are different. The asymmetric cases (a) present, for both the intermediate 
position xq = 10 (solid curve) and the far field position xo = 50 (dashed curve), two 
temporal evolutions. For xo = 10 a local maximum, followed by a minimum, is visible in 
the energy density, then the perturbation is slowly amplifying and the transient can be 
considered extinguished only after hundreds of time scales. For xq = 50 these features 
are less marked. It can be noted that the far field configuration (xq = 50) has a faster 
growth than the intermediate field configuration (x = 10) up to t = 400. Beyond this 
instant the growth related to the intermediate configuration will prevail on that of the 
far field configuration. In the symmetric cases (b) the growths become monotone after 
few time scales (t = 20) and the perturbations quickly reach their asymptotic states 
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(around t = 50). The intermediate field configuration (xq = 10, solid curve) is always 
growing faster than the far field configuration (x = 50, dashed curve). This particular 
case shows a behavior that is generally observed in this analysis, that is, asymmetric 
conditions lead to transient evolutions that last longer than the corresponding symmetric 
ones, and demonstrates that the transient growth for a longitudinal station in the far 
wake can be faster than in the intermediate wake. It should be noted that, even if the 
asymmetric perturbation leads to a much slower transient growth than that observed 
for the symmetric case, the growth rate become equal when the asymptotic states are 
reached (see for example fig. [9]). The temporal window shown in fig. 3 (t = 500) does not 
yet capture the asymptotic state of the asymmetric input. However, we observed that 
further in time the amplification factor G reaches the same order of magnitude of the 
symmetric perturbation. 

The more noticeable results presented in fig. 3 are that the asymmetric growths in 
the early transient are much less rapid than the symmetric ones and that the function G, 
in the case of asymmetric perturbations only, shows a modulation which is very evident 
in the first part of the transient, and which corresponds to a modulation in amplitude 
of the pulsation of the instability wave, see fig|4} In fact, the pulsation varies: in the 
early transient it oscillates around a mean value with a regular period, which is the same 
visible on G, and with an amplitude which is growing until this value jumps to a new 
value around which oscillates in a damped way. This second value is the asymptotic 
constant value. This behavior is always observed in the case of asymmetric longitudinal 
or oblique instability waves. Instead, it is not shown by transversal (0 = vr/2) waves or 
by symmetric waves, see fig. 4, where, on the one hand, the asymptotic value, nearly 
equal to that of the asymmetric perturbation, is rapidly reached after a short monotone 
growth and where, on the other, the growth is many orders of magnitude faster, and as a 
consequence, a modulation would not be easily observable. Thus, we may comment here 
on the fact that two time scales are observed in the transient and long term behavior of 
longitudinal and oblique perturbations: namely, the periodicity associated to the average 
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value of the pulsation in the early transient, clearly visible in the asymmetric case only, 
and the final asymptotic pulsation. The asymptotic value of the pulsation is higher than 
the initial one, typically is about 2.5 times higher. The period of the frequency modulation 
of the energy density G is larger, nearly (1.4 - 1.7) higher, than the average period of the 
oscillating wave in the early transient, because G is a square norm of the system solution. 
Thus the evolution of the system exhibits two periodic patterns at different frequencies: 
the first, of transient nature, and the other of asymptotic nature. When the average 
damping of the energy density in the early transient is not strong and is then followed by 
a monotone asymptotic growth, the change of the frequency of the oscillation is evident, 
see figjij 

This kind of behavior is often observed in the study of linear systems with an oscillat- 
ing norm, a problem that naturally arises in the context of the linearized formulation of 
convection-dominated systems over finite length domains. The occurrence of oscillating 
patterns in the energy evolution of the solutions is linked to the non-normal character of 
the linear operator which describes the system, see e.g. Coppola and De Luca (2006) [30J . 

Figure [5] illustrates an interesting comparison between two-dimensional and three- 
dimensional waves (note that a logarithmic scale on the ordinate is used in part (a) of Fig. 
[5]). The purely two-dimensional wave (solid curve) is rapidly reaching a first maximum of 
amplitude (at about t = 15), then the perturbation decreases while oscillating and reaches 
a minimum around t = 150. Afterwards, the disturbance slowly grows up to t w 300, 
where an inflection point of the amplification factor G occurs. It should be noted that this 
behavior is controlled solely by the evolution of the u z , that is by the vorticity component 
present in the basic flow only. Then, the growth becomes faster and the perturbation is 
highly amplified in time. The purely orthogonal perturbation (dashed curve) is instead 
immediately amplified. The trend is monotone, and does not present visible fluctuations 
in time. The initial growth is actually rapid and an inflection point of the amplification 
factor G can be found around t = 50. Beyond this point, the growth changes its velocity 
and becomes slower, but still destabilizing. Both cases have asymmetric initial conditions 
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and are ultimately amplified. In agreement with Squire theorem, the two-dimensional case 
turns out to be more unstable than the three-dimensional one, as the 2D asymptotically 
established exponential growth is more rapid than the 3D one (see solid and dashed curves 
in Fig. |5](a) for t > 400). However, it should be noted that, for an extended part of the 
transient (up to about t w 380), the three-dimensional perturbation presents a larger 
growth than the two-dimensional one. 

Fig. [6] demonstrates that purely orthogonal three-dimensional unstable perturbations 
may become damped by increasing their wavenumber (k — 7). In the case displayed in 
this figure, this happens when the wave number is increased beyond the value 1. Before 
the asymptotic stable states are reached, these configurations yield maxima of the energy 
density (e.g. when k — 1.5, G ~ 3 at t ~ 30) in the transients. This trend is also typical 
of oblique and longitudinal waves, and it can be considered a universal feature in the 
context of the stability of near parallel shear flows. It should be noted that, in fig. [6j the 
perturbation is symmetric and again the amplitude modulation is not observed in the 
early transient, even in the asymptotically stable situations. However, this example of 
transient behavior also contains a feature which is specific of orthogonal, both amplified 
or damped, and symmetric or asymmetric, perturbations, namely, the fact the most 
amplified component of the vorticity is the u y , see also fig|5j part d. 

In fig. [7] a significant phenomenon is observed for a longitudinal wave. By changing 
the order of magnitude of aj, it can be seen that perturbations that are more rapidly 
damped in space (see, in fig. [TJd, the longitudinal spatial evolution of the wave) yield 
a faster growth in time. In fact, for nearly uniform waves in x direction («j — > 0) 
the configurations are asymptotically damped in time, while for increasing values of the 
spatial damping rate the perturbations are amplified in time (note that a logarithmic scale 
is used on part (a) of the figure). A possible general explanation is that the introduction 
in a physical system of a spatial concentration of kinetic energy is always destabilizing, 
hence the higher is the disturbance concentration the faster is the growth factor. 
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3.2 Physical interpretation of the different growth rate of sym- 
metric and asymmetric disturbances 

The dramatically increased growth rate of the symmetric mode with respect to the asym- 
metric mode can be understood from the induced velocity of the vorticity field. For 
simplicity, imagine that the wake consists of two parallel shear layers of opposite vortic- 
ity, as shown in figure 8(a). Further assume that the vorticity is discretized into a finite 
number of identical vortices. Suppose the upper shear layer is perturbed into a sinusoidal 
shape, as shown in figure 8(b). The induced velocity at the crest of the sinusoid at point 1 
due to the other vortices in the upper shear layer alone is up and to the right as indicated 
by the arrow, corresponding to the classical Kelvin-Helmholtz instability. 

Now consider the asymmetric mode, so that the lower shear layer is a sinusoid in 
phase with the upper one (figure 8(c)). The induced velocity at the crest of the lower 
shear layer at point 2 due to the other vortices in the lower shear layer alone is exactly 
the opposite of that at point 1. Thus the growth in the asymmetric mode is due to only 
higher order gradients in the induced velocity field of one shear layer on the other. 

In contrast, with a symmetric perturbation, the lower shear is the mirror image of the 
upper. In figure 8(d), the induced velocity at point 2 is necessarily the mirror image of 
that at point 1. Both points move in concert to the right, without the low order velocity 
cancellation of the asymmetric mode. Thus the symmetric mode grows much faster. 

This linear absolute instability takes place in the intermediate and far field and acts as 
a source of excitation for the pair of steady recirculating eddies in the lee of the cylinder. 
The onset of a time periodic flow, a supercritical Hopf bifurcation ([31], [32]) indicates 
that in the end the vortices are shed alternatively from the separated streamlines above 
and below the cylinder forming the Von Karman vortex street. The vortex street is 
the stable configuration after the bifurcation has taken place. It has the symmetry of a 
traveling sinuous mode, which indicates asymmetry up-and-down of the cylinder. Thus, 
it can be observed that, also in the context of the vortex shedding, asymmetry shows 
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higher stability properties with respect to symmetry. 

3.3 Asymptotic fate and comparison with modal analysis 

Figure [9] presents a comparison between the initial-value problem and the asymptotic 
theory results represented by the zero order Orr-Sommerfeld problem ([23]) in terms of 
temporal growth rate r and pulsation u. 

In fig. [9]the imaginary part a% of the complex longitudinal wavenumber is fixed, and differ- 
ing polar wavenumbers (k = a r ) are considered. For both the symmetric and asymmetric 
arbitrary disturbances here considered, a good agreement with the stability characteristics 
given by the multiscale near-parallel Orr-Sommerfeld theory can be observed. However, 
it should be noted that the wave number corresponding to the maximum growth factor in 
the case of asymmetric perturbations is about 15% lower than that obtained in the case 
of symmetric perturbations and that obtained by the normal mode analysis. When the 
perturbations are asymmetric, the transient is very long, of the order of hundreds time 
scales. This difference can be due either to the fact that the true asymptote was not yet 
reached, or to the fact that the extent of the numerical errors in the integration of the 
equations is higher than that obtained in the case of symmetric transients, which last 
only a few dozen time scales. Note that this satisfactory agreement is observed by using 
arbitrary initial conditions in terms of elements of the trigonometrical Schauder basis for 
the L 2 space, and not by considering as initial condition the most unstable waves given 
by the Orr-Sommerfeld dispersion relation. Moreover, a maximum of the perturbation 
energy (in terms of r) is found around k = 0.8 and confirmed by both the analyses. 

As shown in fig. [9j we have also contrasted our results with the laboratory experimen- 
tal results obtained in 1989 by Williamson |27j. who gave a quantitative determination 
of the Strouhal number and wavelength of the vortex shedding - oblique and parallel 
modes - of a circular cylinder at low Reynolds number. The comparison is quantitatively 
good, because it shows that a wave number close to the wave number that theoretically 
has the maximum growth rate at Re = 50 (see part a of fig. [9]) has a - theoretically 
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deduced - frequency which is very close to the frequency measured in the laboratory. At 
this point, also the laboratory experimental uncertainty, globally of the order of a ±10% 
in an accurate measurement set up, should be introduced. The uncertainty associated 
to the laboratory method and to the theoretical model (estimated through the difference 
between the position of the maximum growth rate showed by the two cases of asymmet- 
ric and symmetric perturbation) overlaps, which confirms the quality of this comparison. 
The same quantitative agreement is observed also at Re = 100. 

4 Conclusions 

The three-dimensional stability analysis of the intermediate asymptotics of the 2D bluff- 
body viscous growing wake was considered as an initial-value problem. The velocity- 
vorticity formulation was used. The perturbative equations are Laplace-Fourier trans- 
formed in the plane normal to the growing basic flow. The Laplace transform allows for 
the use of a damped perturbation in the streamwise direction as initial condition. In 
this regard, the introduction of the imaginary part of the longitudinal wavenumber (the 
spatial damping rate) was done to explicitly include in the perturbation, which other- 
wise would have been longitudinally homogeneous, a degree of freedom associated to the 
spatial evolution of the system. 

An important point is that the vorticity-velocity formulation, Fourier-Laplace trans- 
formation, allows (over a reasonable lapse of computing time) for the following of the 
temporal evolution over hundred of basic flow time scales and thus to observe very long 
transients. Such a limiting behaviour would not have been so easily reached by means of 
the direct numerical integration of the linearized governing equations of the motion. 

Two main transient scenarios have been observed in the region of the wake where 
the entrainment is present, the region in between xo = 10 diameters (intermediate) and 
xo = 50 diameters (far field), for a Re number equal to 50 and 100. A long transient, where 
an initial growth smoothly levels off and is followed either by an ultimate damping or by a 
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slow amplification for both oblique or 2D waves, and a short transient where the growth or 
the damping is monotone. The most important parameters affecting these configurations 
are the angle of obliquity, the symmetry, the polar wavenumber and the spatial damping 
rate. While the symmetry of the disturbance is remarkably influencing the transient 
behaviour leaving inalterate the asymptotic fate, a variation of the obliquity, of the polar 
wavenumber and of the spatial damping rate can significantly change the early trend 
as well as the final stability configuration. Interesting phenomena are observed. A first 
one is that, in the case of asymmetric longitudinal or oblique perturbations, the system 
exhibits two periodic patterns, a first, of transient nature, and a second one of asymptotic 
nature. A second phenomenon is that, in the case of an orthogonal perturbation, albeit 
always initially set equal to zero, the transversal vorticity component is the vorticity 
component which grows faster. The third phenomenology is linked to the magnitude of 
the spatial damping rate. Perturbations that are more rapidly damped in space lead to 
a larger growth in time. 

For disturbances aligned with the flow, the asymptotic behaviour is shown to be 
in excellent agreement with the zero order results of spatio-temporal multiscale modal 
analyses and with the laboratory determined frequency and wave length of the parallel 
vortex shedding at Re = 50 and 100. It should be noted that the agreement between 
the initial value problem results and the normal mode theory is obtained not using as 
initial condition the most unstable wave given by the Orr-Sommerfeld dispersion relation 
at any section of the wake, but arbitrary initial conditions in terms of elements of the 
trigonometrical Schauder basis for the L 2 space. 
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Figure 3: Effect of the symmetry of the perturbation, (a) - (b): The amplification factor 
G, and (c)-(d) the perturbation vorticity components at y = 1 as a function of time. 
(a)-(c) asymmetric initial condition, (b)-(d) symmetric initial condition. Intermediate 
(x = 10, solid curves) and far field (x = 50, dashed curves) wake configurations. The 
periods T inter ,Tf ar are the periods of the modulation visible on G, in the intermediate 
and far field, respectively. The values of the vorticity component in part (d) have no 
physical meaning. The plot simply shows that, on the contrary of the asymmetric case in 
part(c), the symmetric disturbance growth has a short transient after which it becomes 
homogeneous in time. 



24 



1.5 



0.5 



Re=100 


sym 


k=0.6 


a = 0.02 
i 


asym 


§ = n/4 




x Q =10 







100 200 4 300 400 500 



Figure 4: Pulsation behavior. The wave parameters are those shown in the previous 
figure, for x = 10. 
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Figure 5: Effect of the angle of obliquity 0. (a) The amplification factor G and (b) 
the temporal growth rate r as functions of time. Asymmetric initial condition, = 
(solid curves), = x/2 (dashed curves), (c)-(d) the perturbation vorticity components 
as functions of time at y — 1, (c) = and (d) = ir/2. The periods r inter is the period 
of the modulation visible on G, in the intermediate field. 
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Figure 6: Effect of the polar wavenumber k. (a) The amplification factor G and (b) the 
temporal growth rate r as function of time. Symmetric initial condition, k = 0.5, 1, 1.5, 2. 




Figure 7: Effect of the spatial damping rate (a) The amplification factor G as function 
of time and (b) the wave spatial evolution in the x direction for k = a r = 0.5. Asymmetric 
initial condition, «j = 0,0.01,0.05,0.1. 
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Figure 8: Induced velocity field of the symmetric and asymmetric modes, (a): Idealization 
of the base flow of the wake into two shear layers of opposite sign, (b): Perturbed upper 
shear layer, (c): Induced velocities in the asymmetric mode. Note that they cancel to 
lowest order, (d): Induced velocities in the symmetric mode. Superposition enhances the 
amplitude at both points 1 and 2. 
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Figure 9: (a) Temporal growth rate and (b) pulsation. Comparison among the asymp- 
totic results obtained by the IVP analysis (circles: symmetric perturbation; triangles: 
asymmetric perturbation) and the normal mode analysis (solid lines, see (21] ■ The asymp- 
totic values for the IVP analysis are determined when the condition dr/dt < e (e ~ 10~ 4 ) 
is satisfied. For the symmetric perturbations the asymptote was reached before 50 time 
scales were elapsed, for the asymmetric perturbation the asymptote was determined after 
500 time scales were elapsed. 
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